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1. Introduction 

One purpose of this article is to review the setting for the model of spatial permu- 
tations and its relation with the quantum Bose gas, and to summarize some of the 
material presented in a recent collaboration with Volker Betz 2 . Another purpose 
is to compute the effective interaction between permutation jumps. It involves the 
original interaction potential between quantum particles. While several mathemat- 
ical questions remain unanswered, it is argued that the model of interacting spatial 
permutations describes the quantum interacting Bose gas exactly, and in a simpler 
way. The main phenomenon in bosonic systems is the Bose-Einstein condensation. 
We discuss the links between this phase transition and the occurrence of infinite 
cycles in random permutations. 

Given points xi, . . . ,xn in M'', one considers random permutations tt of iV ele- 



ments with weight 



JV 



4/3 1 

= 1 



Permutation jumps are essentially finite, but permutation cycles can be large. This 
model is illustrated in Fig. [1] It is motivated in large part by the Feynman-Kac 
representation of the Bose gas. We actually discuss a more general setting where 
permutation jumps interact. 

The precise setting is introduced in Section [51 We recall the Feynman-Kac 
representation of the Bose gas in Section [31 it makes the relation between the ideal 
Bose gas and non-interacting spatial permutations clear. The two-body interaction 
between permutation jumps, that is expected to give the exact behaviour to lowest 
order in the strength of the particle interactions, is computed in Section [4l Finally, 
we describe a simple model of interacting permutations in Section [5l It is exactly 
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Figure 1. Illustration for a random set of points x = 
(xi, . . . ,Xn)j and for a permutation tt G Sn- Isolated points are 
sent onto themselves. Permutation jumps are small, but long cycles 
can occur nonetheless. 

solvable, and it provides some understanding about the effects of interactions on 
the Bose-Einstein condensation. 



2. The model of spatial permutations 

Let A C M'' be a cube of size L and volume V = L''^, and let N eN. The state 
space of the model of spatial permutations is 

Qa^n = A"^ X Sn, (2.1) 

with <Sjv the symmetric group of permutations of N elements. We are interested 
in the properties of permutations, and all our random variables are functions 9 : 
Sn — > Their probability distributions depend on spatial variables in an indirect 
but essential way. Let iiin) denote the length of the cycle that contains i, i.e. the 

smallest integer n > 1 such that n"'{i) — i. The most important random variable 
is the density of points in cycles of certain lengths. For n, n' e N, let 



1 



e„,„' W = yif^ii = l,...,N ■.n< tiiir) < n'). 

The expectation of the random variable 6 is defined by 

1 



(2.2) 



Sa,jv(^) 



Z{K,N)N\ 



f dxY e{Tz)e-"^'^'^K (2.3) 



Here, the normalization factor Z{K,N) is chosen so that E\^n{^) = 1- The term 
N\ is present in order that Z{K, N) scales like the exponential of the volume of A 
— thus behaving like a partition function in statistical mechanics. The integral is 
over N points in A, denoted x = (x\, . . . , Xn). 
We consider Hamiltonians of the form 



N 



H{x,Tt) =^i{x^ - X^^i-j) + ^ V{Xi,X^(^i-),Xj,X^y^), 



(2.4) 



i=l 



l<i<j<N 
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with ^ a spherically symmetric function R"* — )■ [0, oo], and V a translation invariant 
function R"*** — > R. We also suppose that £, is increasing and that ^(0) = 0. One 
should think of typical permutations as involving finite jumps, i.e. \xi — ^^^(i) | stays 
bounded as A, — > oo. 

The major question concerns the occurrence of infinite cycles. It turns out that 
the distribution of cycles can be well characterized in the absence of interactions, 
with the potential V = 0. We need a few hypotheses on ^. Let C = J . 
We suppose that e~'> has positive Fourier transform, which we denote C e^^^'^^ . 
Precisely, we have 

Ce-'^''^ = / 6-2^^"=^ e-«(^) dx. (2.5) 



The case of physical relevance is £,{x) — jplx]"^ with j3 the inverse temperature, in 
which case e{k) — 47r^/3|/c|^. But it may be of mathematical interest to consider 
other functions, including some where has bounded support. Criteria that 
guarantee positivity of the Fourier transform are discussed e.g. in [5]. 
We define the critical density by 

The critical density is finite for d > 3, but it can be infinite for d ^ 1,2. The 
experienced physicist will have recognized the formula for the critical density of 
Bose-Einstein condensation. The relation with the Bose gas will be discussed in 
the next section. In the following theorem we fix the density p and we let N = pV 
in the expectation (|2.3p . 



THEOREM 1. Let ^ satisfy the assumptions above. Then for any Q < a < h < 1, 

and any s > 0, 

(a) hm £^A,pv(£>i,ya) = < 

(b) lim Ea,pv{Qv- yO = 0; 

{0 ifp<Pc; 
s ifO<s<p-pc; 
p- Pc ifO<p-Pc<s. 

In order to understand the meaning of these claims, one should think of a as 
barely bigger than 0, and b barely smaller than 1. In part (a), Qi ya is the density of 
points in finite cycles. All points belong to finite cycles if p < Pc- However, if p > pc, 
a fraction p — pc of points belong to infinite cycles. It is natural to ask oneselves 
about the size of "infinite cycles" in a finite domain of volume L'^ = V . One could 
expect the typical length to be of order L^, since the continuum limit of random 
walks has Hausdorff dimension 2, and cycles are somewhat like closed random walks. 
However, part (b) shows that cycles of length , with < 6 < 1, have vanishing 
density. Thus infinite cycles are macroscopic, i.e. each cycle involves a strictly 
positive fraction of points. The statistics of macroscopic cycles is characterized in 
part (c). 

The proof of Theorem [1] can be found in ^ . Actually, the correct statement 
involves periodic boundary conditions; the interested reader is invited to look in [2] 
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for the precise statement. Theorem [T] extends an earher result of Siito for the ideal 
Bose gas [13]. 

A different model is investigated in [1], where the positions Xi,.. .,xn form a 
cubic lattice. The density is always equal to 1, but the function ^{x) depends on a 
parameter that represents the temperature of the system. It is found numerically 
that the critical temperature for the occurrence of infinite cycles is close but different 
from that of the ideal Bose gas. Many properties are similar, however; infinite cycles 
are also macroscopic. A surprising fact is that the expectation of the length of the 
longest cycle seems to be identical to that in the ideal Bose gas; this suggests that 
the distribution of macroscopic cycles may be the same. 



3. Feynman-Kac representation of the Bose gas 

The Feynman-Kac formula relates the kernel of e~'^^ , with H a Schrodinger 
operator, to the Brownian motion, whose mathematical expression is the Wiener 
measure. It seems to have first appeared in [3], precisely in the context of bosonic 
systems and in the discussion of cycles. Ginibre wrote an excellent mathematical 
introduction to the Wiener measure, the Feynman-Kac formula, and its application 
to bosonic systems [S]. We review these notions here without introducing the full 
mathematical setting, but all equations below can be justified with a bit of analysis. 
In particular, we do not discuss the details arising from the boundary conditions; 
as usual in statistical mechanics, they are irrelevant for large systems. 

Let A C M'' be a cube of size L, and let denote the normalized Gaussian 
function 

It is not hard to check that 

gsix - a)gt{x - b)Ax = gs+t{a - b), (3.2) 

a 

and that, after iteration, 

P n 

/ dX2 ■ ■ -dXriY^gtiXi+l - Xt - Gi) ^ gnt{xn+l - Xi - Y^'l^j^Qt). (3.3) 

Let H = — A + t/bea Schrodinger operator in L^(A), with A the Laplacian and 
U a smooth real function. This operator is unbounded and we need to specify its 
domain. We can choose the space of functions on A with Dirichlet boundary 
conditions. Then H is symmetric and we consider its self- adjoint extension. Of 
relevance to statistical mechanics is the operator e~^^ . It is a nice operator, 
bounded and compact, but these properties are not important here. The Feynman- 
Kac formula states that (with Xn+i = xi) 

Tr e"''^ = lim / dxi . . .dxn\Y\_92i3/n{xi+i - Xi) expj y^^U{xi) \ 

dxexp{-i^ U{^{s))dsY'^fA^)- 

Here, w is a Brownian bridge starting and ending at x and traveling in time 2/3, and 
Wff is the Wiener measure. In the second line we should restrict the paths to stay 
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inside A, because of Dirichlet boundary conditions; we neglect these technicalities, 
however. 

Let us turn to the description of bosonic systems. The state space for N quan- 
tum bosons in a domain A C M'' is the subspace Lgyj-^(A^) of symmetric complex 
functions of variables. The Hamiltonian is given by the Schrodinger operator 

N 

ir = -^A, + J2 U{x,-xj). (3.5) 

i=l l<iJ<N 

Here, Ai denotes the d-dimensional Laplacian for the i-th variable, and U{xi — Xj) 
is a multiplication operator that represents the interaction between particles i and 
j. We always suppose that U{x) > 0. We can choose the self-adjoint extension 
of H that corresponds to Dirichlet boundary conditions. Of course, the sum of 
A^ Laplacians in A'' can be viewed as a Laplacian in A''^, so we can apply the 
Feynman-Kac formula. 

The canonical partition function is equal to 

TvL2^^^A-)e-^" =Tri2(A«)P+c-''^ (3.6) 

where P-f is the projector onto symmetric functions, 

P+if{xi,...,XN)^ ^ (/?(x^(i), . . . ,a;^(Ar)). (3.7) 

tt^Sn 

Using this projection and the Feynman-Kac formula, the partition function of 
the Bose gas can be written as 

(3.8) 

with the Gibbs factor given by 

N „ „2/3 

^ J-Q / ^wZ^^^^iu^^)^ exp{-i / U{u:.[s)-u;,{s))ds]. 

1=1"^ l<i<j<A''^'' 

(3.9) 

This formula is illustrated in Fig. [2l It involves spatial positions and permutations 
of these positions. 

In the case of the ideal gas, U = Q, the paths no longer interact and the Wiener 
integrals in (|3.9[) can be computed. We find that 

with 



i/W(a.,7r)^i-f]|x.-x^(,)|'. (3.11) 



' 1=1 

The prefactor in p.lOp plays no role in expectations of random variables and it can 
be ignored. Thus the ideal Bose gas is equivalent to the "ideal" model of spatial 
permutations with Gaussian weights. Random variables of permutations have same 
distribution in both models, and the transition to infinite cycles takes place at the 
same critical density. 

The equivalence between the occurrence of infinite cycles and Bose- Einstein con- 
densation is an open problem. It is known to be true in the ideal gas, see p ^ fTS l fTS ] . 
but it does not seem to be true in strongly interacting systems in a solid phase. 



DANIEL UELTSCHI 

P 




Figure 2 . Feynman-Kac representation of a gas of N bosons. The 
horizontal plane represents the d spatial dimensions, and the ver- 
tical axis is the imaginary time dimension. The picture shows five 
particles and two cycles, of respective length 4 and 1. 



Pollock and Ceperley have argued that superfluidity is related to spatially winding 
cycles [Ih - Such cycles are clearly infinite in the thermodynamic limit. On the 
other hand, we know from Theorem [T] that infinite cycles are macroscopic (i.e. they 
have strictly positive density), so they certainly have non-zero winding number. 
Infinite and winding cycles should therefore be equivalent. Superfluidity is by no 
means equivalent to Bose-Einstein condensation. These facts bring some level of 
confusion and we can only hope that they will be clarified in the future. 

However, it is expected that, in dimension c? > 3, weakly interacting bosonic 
systems have the same critical density for Bose-Einstein condensation, superfluidity, 
infinite cycles, and winding cycles. Hereafter, we study the occurrence of infinite 
cycles in the weakly interacting regime, and we implicitly assume that they reveal 
a Bose-Einstein condensation. 



4. Exact two-body interaction for permutation jumps 

4.1. Expansion of path interactions. The two-body interactions between quan- 
tum particles translate into many-body interactions for permutations. But we can 
perform an expansion and see that, to lowest order, we obtain a two-body interac- 
tion between permutation jumps. 

LetW^y — g^^{x—y)W^y be a Wiener measure normalized such that / dW*_y{uj) 
1. From (13.91). we have 



with given by (I3.1ip . and 

TV 



1=1'' l<i<j<N 
N „ 

= [n/*i..Mn(i-T^M)- 
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The last product is over bonds b — with i ^ j. For i < j, we defined 

ujb = LOi — LUj , and 

T{LUb) = 1 - c-iJ'o" UMs))ds (4 3) 

Expanding the product in (|4.2p . we have 

N ^N(N-l) k 

e-^'"(-^' = [n/dM?.'i.,M E (-1)' E n^^KJ- (4-4) 

j=l fc=0 {bi,...,bfc} r?x=l 

In the regime of weak interactions, the typical k in the above sum is a small fraction 
of the volume, and the typical 5i, . . . , 6^ are mostly disjoint. 

We first perform the expansion in a somewhat cavalier fashion. We will be more 
precise in Section where we will check that we have identified the leading order. 

Let 6„i {irrnjm}- Let us assume that bi n = for all ^ 7^ m; then 

JV „ fc fc 

[n / n ^(^fc™) = n yi^^^^^^i^^j^^j^^^^u,.)) (4-5) 

z— 1 m—1 m—1 

where the potential V has been defined by 

V{x,y,x',y') = [ dW^^iu;) [ d1?,%,(c.')T(c. - c'). (4.6) 



This is the two-body interaction between jumps x >—>■ y and x' ^ y' . The expression 
(|4.6p can be simplified, see Eq. (|4.24p below. We use the identity (|4.5p for aU 
bl,. . . ,bk that appear in (j4.4p . not only disjoint ones. This is an approximation; 
it assumes that either the terms with intersecting b^'s are not important, or that 
their contribution is close to (14.51). We obtain 



k 

g-if<^)(.,.) ^^(_l)fc ^ J]y(a:,„,x.(,„),x,„,a;.o„)). (4.7) 

fc>0 {bi,...,bfc} m=l 

Ignoring the possibility that a same bond may occur several times, we get 

fc>0 ' bi,...,bkm=l (4-8) 

= exp|- ^ V{xi,x^(^,),Xj,x^(^j))y 

l<i<j<N 

These approximations suggest that, to lowest order in the strength of the interac- 
tion, the multi-body interaction arising from the Feynman-Kac representation can 
be approximated by the two-body interaction defined in (|4.6p . 

4.2. Cluster expansion. It is not clear that the approximations above have pro- 
duced the correct terms, that are exact to lowest order. In this section we perform 
a cluster expansion. It cannot be entirely justified from a mathematical point of 
view, but it nevertheless clarifies the approximations. 

Consider the graph with vertices {1, . . . , /c}, and with an edge between £ and m 
whenever b^nbm 0- We say that a set of bonds B = {6i, . . . , 6^} is connected ii this 
graph is connected. Let suppS = Ubgsb. We say that B and B' are compatible 
if their supports are disjoint, suppB fl suppB' — 0. Then the sum over sets of 
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bonds in ()4.4|) can be written as a sum over connected and mutually compatible 
B's, namely 

^■>0 m^l £>0 ■ Bi,...,Bi m^l 

compatible 

(4.9) 

The contribution of compatible B's factorizes. For a connected B, let us introduce 
^Bix,n) = (-1)1^1 [ n / dW^'l^l^O] n (4-10) 

Notice that ^b{x,7t) depends only on positions Xi and Xjrfi) for i g suppB. Then 
we have 

y: ti'i'B^M. (4.11) 

^>0 ■ Bi,...,Bi m=l 
compatible 

We now apply the cluster expansion method, see e.g. [SI El] for references. Given 
Bi, . . . , Bi, let ip{Bi, . . . , Bi) be the following combinatorial function: 



1 if^=l; 



Gce(Si,...,s 



Here, G{Bi, . . . , i?^) denotes the graph with £ vertices, and with an edge between i 
and j whenever Bi and Bj are not compatible. The sum is over all connected sub- 
graphs of i vertices, and \G\ is the number of edges of G. Notice that <f{Bi, . . . , Bi) 
is zero unless Bi, . . . ,Bg form a cluster, i.e. unless G{Bi, . . . , Bi) is connected. 

The cluster expansion yields a convergent series for the logarithm of (|4.1ip , hence 
fori?*^^^. Precisely, 

e 

i/«(a;,7r)--^ ^ (^(Si, . . . , S,) [] (a;, tt). (4.13) 

£>1 m=l 

Let be distinct indices. The previous equation suggests to define the 

fc-body interaction by 

m 

yW((a;,,,x,(,,))tO =-E E V{B^,....B^)\{^B,M- 

m>l Bi,...,Bm 1=1 

Ui supp Bf = {ii ,...,ifc} 

(4.14) 

Then i/^^^ is given by 

H('Hx,^) = Y. E ^^'nK.^-(^.))ti)- (4-15) 

fc>2 l<!:i< - <ifc<Af 

Everything here is exact, and it is rigorous provided we can prove the absolute 
convergence of the series of cluster terms in (|4.13p . A sufficient criterion is that, 
for any i, 

|$B(a;,7r)|e''l^l <a (4.16) 
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for some constant a > 0. See e.g. [9l[T4] for concise statements about cluster expan- 
sions. The sum above involves bonds whose positions are far away. In order to get 
such an estimate, one needs to control spatial decay. It depends on permutations, 
and there are combinatorial difhculties. 

We conclude this subsection by discussing various estimates for the terms above. 
Using 1 — e~^ < x, we have that 

||T|U = supT(w) </3||C/|U. (4.17) 

The interesting regime of parameters is P and J7 ^ 0, so T is arbitrarily 

small. If the potential [/ is a hard-core with small radius, then ||T||oc — 1, but 
||T||p is small for p < oo. We also have that, for any B, x, and tt, 

|$B(a;,^)|<||T||^^l. (4.18) 

Consider the series (|4.14|) for the potential at order k. The sets . . . , Bm that 
contribute to lowest order are such that \Bi\ = k — 1. It follows that V^''^ is of 
order ||T||^-i. 

We can extract the lowest order term. The expression for V^^'^ {xi , x^^^i) , Xj , x^^^j) ) 
involves terms of arbitrary orders. But we only need to consider — $3 with B 
containing the single bond h = We then obtain the potential V defined in 

(IMI). 



4.3. A simpler expression for the interaction. We now seek to simplify the 
formula (j4.6p . Namely, we can replace the two integrals over Brownian bridges by 
a single integral, which will lead to the nicer formula (|4.24l) . We have 



dW^^{Lo) / dwX' (w')T(w - to') = lim / dx2 . . . dx„ dx'^ ...dx'^ 

\J\92^3/n{x^+l- Xi)g2fj/n{x[j^i- x'i) 1 - exp| - - ^ C/(a;, - x ■) | (4.19) 



n 

i=l 



with xi ~ X, Xn+i = y, x[ = x' , x'^^^ — y' . Let us introduce Zi = Xi — x[. It is not 
hard to check that 

g2t3/n{xi+i - Xi) g20/nix'i+i - x'^) = - Xi - ^Z^+i + \zi) g4(3/n{zi+l " ^i) . 

(4.20) 

Substituting into (|4.19p . we get 

J ^ ■ ■ •dz„|^J|540/„(zi+i - Zi)j 1^1 - exp|--- ^ J7(zi)| 

» n 
/ dX2 ■ ■ -dXriY^gfl/niXi+l - - ^Zi+i + ^Zi). (4.21) 

Using Eq. I\3.'6\i . the last line is equal to 

90{y-x - ^{y - y' - X + x')) ^ gp{^{y + y' - x - x')) . 

The first line of (|4.2ip yields an integral over Brownian paths. Then (|4.2ip is equal 
to 

g^{liy + y'-a:-x')) [ [l " e"^ ^ ^(-(^»<^^ ] dW^f (c.). (4.22) 
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Finally, we have the following identity, similar to (I4.20p 

52/3(y - x)g2p{v' - x') = 5/3(^(2/ + y' -X- x'))gip{y - y' - x ^ x'). (4.23) 

Recall that the two-body interaction defined in (|4.6|) involves normalized Wiener 
measures. Putting normalizations back, we get the following elegant formula for 
the interaction between jumps x ^ y and x' ^ y' , 



Vix, y, x', y') = 111- c-i /o" ^(-(^))^1dW^,^^,, ,_,,(.). 



(4.24) 



It would be useful to obtain a closed form expression in terms of special functions, 
if it is possible. When U consists of a hard-core potential of radius a, V{x, y, x' , y') 
is equal to the probability that a Brownian bridge, starting at x — x' and ending at 
y — y' , intersects the ball of radius a centered at 0. 

4.4. EfTect of interactions on the critical temperature. The model of spa- 
tial permutations should help to clarify the effects of interactions on the critical 
temperature of Bose-Einstein condensation. 

Let Ti"^ be the critical temperature for Bose-Einstein condensation as a function 
of the scattering length a of the interaction potential U between quantum particles. 
It is believed that Ti"^ behaves in three dimensions as 



-ho(pi/3a), (4.25) 



with c a universal constant that does not depend on the mass of particles or on the 
interactions. The value and even the sign of c has been contested in the physics 
literature, although a consensus has recently emerged that c « 1.3. See [II [71 [8l flO] 
and references therein. 

The model of spatial permutations is clearly simpler than the Feynman-Kac 
representation of the Bose gas, and is therefore better suited to Monte-Carlo simu- 
lations. More importantly, we expect that this model, with the interaction (|4.24p . is 
exactly related to the original quantum boson model, to lowest order in a. Numerical 
simulations should allow to determine the constant c in the model of permutations 
with high precision, and with high confidence. It should be identical to the universal 
constant of (I4.25|) for the interacting Bose gas. 

5. A SIMPLE MODEL OF INTERACTING SPATIAL PERMUTATIONS 

In this final section, we discuss a simple model of interacting spatial permuta- 
tions that was introduced in [2J . We consider only interactions between permutation 
jumps of 2-cycles, arguably the most important. The resulting model is exactly solv- 
able, and it provides a heuristic description for the shift in the critical temperature 
of the Bose-Einstein condensation. 

The approximation consists in replacing the Hamiltonian ()2.4|) by 

1 ^ 

H{x,Tr) ^ —^\x^ - X^l^.^l'^ + ^ V{X^,XJ,XJ,X^). (5.1) 

^ 1=1 l<i<j<N 

The interaction term V{-) is given by (|4.24p as before. From now on, all com- 
putations will be exact, at least to lowest order in the scattering length of the 
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original potential U . We consider the three-dimensional case, obviously the most 
interesting. A computation shows that 

V{x,v,y,x)^-^ + 0{a'). (5.2) 

F - y\ 

The lowest order term in the right side does not depend on /9, surprisingly. The 
expectation of a random variable of permutations is given by (j2.3p , 

i^A,.W-^(^E^wX„d.e-(-). (5.3) 

We now substitute H with the following simpler Hamiltonian i?^"): 

1 ^ 

H(^\x,t:) ^ Y.\x,- + aN2{^), (5.4) 

i—l 

with N2{'k) denoting the number of 2-cycles in the permutation tt. The substitution 
is exact provided that, for any given permutation tt, 

dxe'"^''^''^ = [ da;c--f^'°'(=^^'^) . (5.5) 
Ja" 

Isolating the contribution of 2-cycles, this equation reduces to 
With V{-) in jSS]), we find that 



A2 JA2 



^^)'''« + 0(a^). (5.7) 

We now compute the pressure of the model with Hamiltonian H'^°'\ The grand- 
canonical partition function is given by 

^'(/3,A,M)=E^ E J-^'" (5.8) 

N>0 ■ ireSjv^" 

It is convenient to work in the Fourier space. Let us introduce a new partition 
function, 

Z{P,A,^^) = E ^ E E ^""""^^^^ n (5.9) 

Ar>0 ' fci,...,fciv6A* ir65iv i=l 

Here, A* ~ jT? is the dual lattice. The thermodynamic pressure is defined by 

p("H/?,A^) = Jim -^logZ(/3,A,/i) (5.10) 

One can verify that the partition functions Z and Z' differ in two respects only. 
First, a normalization is missing, which results in a shift of the chemical potential. 
Second, Z has been defined with periodic boundary conditions, unlike Z' (where 
boundary conditions are neither periodic, nor Dirichlet). But both partition func- 
tions yield the same thermodynamics; precisely, 

^ lim log Z'(/3, A, ii) = (/3, A/ + I3 log(47r/3)) . (5.11) 

V— >-oo pV ^ 
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We now compute p^"-'. Introducing occupation numbers, (15. 9|) becomes 



z(/3,A,Ai)= n 



-/3(|27rfc|^-A')nfc 



— QAf2(7rfc) 



(5.12) 



We decomposed the permutation n into permutations (ttj,) for each Fourier mode, 
and we also used 

- E ^2(vrfe). (5.13) 

fceA' 

Notice that the chemical potential needs to be strictly negative, as in the ideal gas. 
We get 



p(")(/3,M) = ^lhn^^El°4E 



-/3(|2^fc|2-M)» 



fceA* '-n>0 



(5.14) 



Let us compute the bracket above. For given tt G iS„, let rj denote the number of 
cycles of length j. Then jfj = n, and the number of permutations for given 
(rj) is equal to 



The bracket in (|5.14|) is equal to 



El x ^ n\ 

n>0 ' ri,r2. - >0 '■'■3>^ ■' 

iZi jrj=n 



' 3- 



Y J_|-l g-i/3(|277fc| = -A.) j'-j ^ J_|-l ^-2/3(|27rfcp-;.)-a]'^2 



ri,r3,r4,--->0j = l,3,4 

exp 



1-2 >o 



exp 



I ^ 1 g-j/3(|2^fcf"-p) _^ 1 g-2/3(|27rfc|"-M)-a I 
j=l,3,4,... 

{- log(l - e-'5(|2^'=l'-^) ) - i e-^^d^-'^l'-^) (1 - e-" )}. 



This can be inserted into (|5.14p . In the limit ^ 00 the expression converges to 
a Riemann integral. If a = 0, we get the pressure of the ideal gas 



P 



(/3,M) = -^ log(l-e-^(|2-'=l^-'^))dfc, 



as expected. And if a 7^ 0, we get 

p(")(/3,/i)=p(°^(/3,M) 



211/27^3/2^5/2 



(1-e— 



(5.15) 



(5.16) 



The pressure p'-"-' is plotted in Fig. [3] (a) as a function of /i. One can consider 
other thermodynamic potentials as well. Recall that the free energy Z*^") is function 
of the (inverse) temperature and of the density, and it is related to the pressure by 
a Legendre transform: 



/(")(/3,p)=sup[pA^-p(")(/3,M)]. 



(5.17) 
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(a) 



(b) 



Figure 3. The pressure and the free energy of the simple inter- 
acting model in three dimensions. 

One then obtains the graph depicted in Fig. [3] (b) . It is strictly decreasing up to the 
critical density pi"'' — pf'\(3), and it is constant afterwards. The critical density is 
equal to the derivative of p^"^ with respect to /i at 0—. We have 

1 



(a) 



(0) 



29/2^3/2^3/; 



■(1 



(5.18) 



The first term of the right side, is equal to the critical density of the ideal gas, 
Eq. (|2.6p . The second term is the correction due to our simple interaction. 

We see that pc"^ is smaller than p"c^ — interactions favour Bose-Einstein con- 
densation. This observation is in line with physicists' expectations. The heuristics 
is particularly simple in this model: 2-cycles are penalized and this favours all other 
cycles, including infinite cycles. The latter occur therefore at a lower density. While 
elementary, this heuristics may well be correct. 

Let us now estimate the change in the critical temperature. Using (|5.7p with 
we find that, to lowest order, 



(a) 



(0) 



(0) 



3^2 



(5.19) 



with c = 0.37. This formula can be compared to (|4.25[) . If we believe the value 
c = 1.3 found numerically, then 2-cycle interactions account only for a fraction 
of the effect of all interactions. One could also take into account the interactions 
within 3-cycles and longer cycles; the constant c would increase a bit. 

One would expect infinite cycles to occur for all densities larger than the critical 
density (|5.18p . More precisely. Theorem [1] should remain valid for a > 0, replacing 
Pc by pc"^. But only a weaker claim has been proved so far. 

THEOREM 2. For anyb<l, 

4 



lim EA^pviQvKpv) > P 



(1 



Theorem [2] guarantees the existence of macroscopic cycles for large enough den- 
sities. The proof can be found in f^. 
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